Magnetoplasmons excitations in Graphene for filling factors z/ < 6 
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In the frame of the Hartree-Fock approximation, the dispersion of magnetoplasmons in Graphene 
is derived for all types of transitions for filling factors v <6. The optical conductivity components 
of the magnetoplasmon curves are calculated. It is shown that the electron-electron interactions 
lead to a strong re-normalization of the apparent Fermi velocity of Graphene which is different for 
different types of transitions. 
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I. INTRODUCTION 

Graphene is a monolayer of Graphite with a band 
structure composed of two cones located at two in- 
equivalent corners K and K' of the Brillouin zone at 
which conduction and valence bands merge. This com- 
pound has recently received a large attention because 
of the unusual sequence of quantum Hall states it re- 
veals [H, 0] ■ In contrast to conventional two-dimensional 
electron gas (C2DEG) which display a quadratic dis- 
persion law, Graphene exhibits a linear dispersion law 
-^(~P*) = ivi?|~p*| as a function of the momentum ~p 
leading to a Dirac's type Hamiltonian with a Fermi veloc- 
ity vf replacing that of the light. Different one-electron 
band structure models, not including electron-electron 
interactions, lead to values vp — 0.86 x lO^m/s with some 
variance, but this is the value which will adopted in this 
report. This peculiar dispersion law has two important 
consequences in contrast to C2DEG (see for instance Q ) : 

(i) the wave functions have a spinor type character and 

(ii) under a magnetic field B applied perpendicular to 
the Graphene plane, the Dirac energy spectrum evolves 
into Landau Levels (LL) with energies given by: 



En — sgn{n)v F \/2ehB\n\ = sgn{n)Eio 



(1) 



where n scans all positive and negative integer values 
including zero. 

Magnetoplasmons (MP), in a two-dimensional electron 
gas, are excitations between LL, known to be described 
in terms of excitonic transitions due to electron-electron 
interactions (EEI): they reveal a specific dispersion as a 
function of the two dimensional wave vector k of the ex- 
citon. For a C2DEG, the theory, derived in the frame 
of the Hartree-Fock (HF) approximation, has been first 
developed [1, 0] for integer values of the filling factor 
V = Ns^o/B {Ns and $o being the two-dimensional car- 
rier concentration and the flux quantum respectively). 
These studies have been extended to the case of non- 
integer values of ^ Q and have also included the cal- 
culation of matrix elements for the optical conductivity 



The effects of EEI in Graphene have recently been 
reported on a theoretical basis @| but with a different 
model than that of Refs. [1, [1, [y| and restricted to in- 
teger values of the filling factor. We have followed here 
the lines of Ref . , which has been shown to reproduce 
quantitatively experimental results foj when they are 
interpreted in terms of MP excitations. 

Because of the Kohn's theorem [l^, the EEI effects 
turn out to be tiny for C2DEG. However this theorem 
does not apply for a linear dispersion law and there- 
fore EEI are expected to induce significant effects in 
Graphene. Indeed recent experimental inve stig ations of 
the magneto-optical transitions in Graphene [llf, have 
been interpreted with an effective velocity c, replacing vp 
in Eq.l, ranging between 1.03 to 1.18 x 10®m/s and show- 
ing a re- normalization of vp which here, we will show, is 
mainly due to electron-electron interactions. 

On general grounds, the Coulomb energy character- 
istic of electron-electron interaction in magnetic field is 
Ec — c^/kIb where k is the electronic dielectric constant 
of the material and Ib = {eB)~^/'^ the magnetic field 
length. The magnetoplasmon approach assumes that Ec 
is smaller than the one-electron energy transitions. In 
the present case of Graphene Ec{meV) = ll.2y/B{T) 
and Eio{7neV) = "H^/WT) (see Eq.l) leading to the 
ratio Eio/Ec = 2.77 that is a condition better fulfilled 
for Graphene than for GaAs based C2DEG and fur- 
thermore not dependent on the value of the magnetic 
field. Magneto-excitons should be therefore more stable 
in Graphene than in C2DEG and the approach derived 
for these later compounds should be valid. The report 
is organized as follows: we will first describe the general 
formalism used to derive the MP dispersion curves (sec- 
tion II). We apply it to the case of filling factors i/ < 2 
in section III and to the case of 2 < ^ < 6 in section IV. 
Results will be discussed and compared to experimen- 
tal results in section V. The details of calculations are 
reported in the appendices. 



II. 



GENERAL FORMALISM 
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In contrast to the GaAs case, Graphene has two valleys 
which lead to the conclusion that in the absence of spin 



2 



splitting and valley splitting each Landau Level (LL) is, 
in general, four times degenerate. The fourfold degener- 
acy of the n = LL is still due to spin and valley symme- 
tries, but two of these levels have an electron-like charac- 
ter and the two other ones a hole-like character. We will 
restrict our analysis to the zero temperature case. Be- 
cause of the peculiar symmetry of the problem the wave 
functions have a spinor character which can be expressed 
in the Landau gauge, with the potential vector compo- 
nents of the magnetic field = ~ and Ay ~ Bx, 
as: 



-isgn{n)ip\n\-i{x - p) 



p) 
p) 



(2) 



where ~p is the two-dimensional vector of components x 
and y, c„ — 1 for n — and l/\/2 otherwise whereas 
sign(n) = 1,0,-1 when n > 0,= 0, < respectively. 
Lp\n\ (x) is the standard normalized Landau wave func- 
tions. Note that these wave functions differ from those 
proposed in Ref. 3 by the phase factor due to the dif- 
ferent gauge used. Following the lines of Ref. ^ , we call 
-^n n' a i ^^"^ Creation operator of an exciton of energy 

Eexik), corresponding to a transition from LL n' with 
spin a in valley i to a LL n of the same spin and same 
valley. This operator is defined as a function of and 
a\, the standard one particle creation and annihilation 
operators respectively, as: 

Al„' . ^(fc)|0) - E exp(ifc.(p + ky/2)) 

' ' P (3) 

The total Hamiltonian of the system is written as: 

Htot = E ^'m<^m,p,a,i^m.p,aA + Hint (4) 
m.p.(T,i 

where fiwj^* is the corresponding one electron energy 
of the LL m with spin a in valley i. The Coulomb inter- 
actions appear in Hint as : 



Hint = \ j dpidp2V{pi - P2) 

where denotes the scalar product and 

Fa4~P) =^F^p{~p)an,p,a,i 



(5) 



n,p 



V{pt -p^) = 



V{q) 



(6) 



with V{q) being the 2D Fourier transform of the 
Coulomb potential V{r) — /[kt). 

After some calculations we obtain an analytic expres- 
sion for Hint which is expressed as: 

Hint = ^ E exp(ig^(pi - P2 - qy)) X 

Jn4,ni{q)Jn3,n2{ ^ 



(7) 



In Eq.7, the summation is extended over the en- 
semble ni, 77.2, 713,714 of LL, the ensemble pi,P2 of the 
y-component of the momentum, the ensemble of spin 
CTi, (72, both valleys ii and «2 and the wavevector q. 

The function J,„ „ is defined as: 



CnCm{sgn{m)sgn{n)Ji„^\_i^\n\-i{q) 



|m|,|n| 



(8) 



with the usual definition of the integral Jm.n{q) valid 
for m > n: 



J m.7i {q) 



iqx 



^ 2 ' 



(9) 



where X™~"(a;) are the Laguerre polynomials. For 
m < n the relation Jm.n(q) = Jn,mi~'l) holds. 

Using the random phase approximation (RPA) to treat 
the combination of creation and annihilation operators 
we arrive to the following expression for the Exciton en- 
ergies (the notation |0) representing the ground state of 
the system) ; 



n' ,n2 ,n' ,712 

(0)- 

(0)]/„^2,. 

xA+ , .(fc)|0) 

{ky^ kx){fnS ~ /n',i)^n2,n4,o-, 7(^)1*^) 



"2 

E^^ 



n' ,712 ,71,714 ' 



E 



^7i' ,7i2 ,''^3 /'I ( ^y) 



27r 



\-i n.% J 71 ,7/ 



(10) 



where j is the filling factor of LL n with spin a in 
valley i. The matrix elements E and V are given by: 



V 



"1 ,"2, "3, "4 



{q) = V{q)Jni,m{q)Jn^,nA-q) 



Fni,n2,n3,ni(,k) 



dq 



■Vn„n2,n,,nAq)e"'-' (H) 
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We note, at that level, that Eq. 10 is formally equiva- 
lent to that obtained for C2DEG @ except for the defi- 
nition of the different matrix elements which here takes 
into account the spinor character of the wave functions. 
In this equation, the second term (second line) is a mea- 
sure of the difference of exchange energies of the LL n' 
and n. The third line is related to the direct electron-hole 
Coulomb interaction (Exciton binding energy). Both 
terms are involving excitons of same spin and same val- 
ley. The last term of Eq. 10 describes the simultaneous 
annihilation and creation of excitons at different points 
of the Brillouin zone (RPA contribution): it includes all 
possible transitions without restriction to spin or valley 
indices. The exchange terms deserve a special attention 
in the present case. The corresponding expression for the 
exchange in Eq. 10 reads as: 



En^m,n,m{0) = |Cnp|Cmp\/2 J dxc 

{L\m\)L\n\ + {sgn{n)sgn{m)f L\rn\-iL\n\ 
2x2 



sgn{n)sgn{m) 



-^l^ml-l-^H-l} (12) 



where all Laguerre polynomials have arguments . In 
Eq.lO (second line), the summation over n2, for these 
terms, has to include all LL from — oo to +oo. The eval- 
uation of the exchange contributions to the different sit- 
uations is given in appendix A. 

Solving Eq.lO results in diagonalizing an Hamiltonian, 
the size of which depends on the number of transitions 
which are assumed to be coupled by electron-electron 
interactions. In reality this number is very large for 
Graphene due to the existence of interband transitions 
but since Ec is smaller than the energy of transitions, we 
can reasonably assume, in the spirit of the HE approxi- 
mation, that EEI will not couple, at first order, transi- 
tions with different one-electron energies. In that case, 
the problem reduces to solve the Hamiltonian for each 
type of optical transitions which depend on the value of 
the filling factor. However, even if one solves the whole 
problem in successive steps, one has to keep in mind that 
a common energy scale should be adopted for all transi- 
tions in order to compare such results with experimental 
ones. 

When writing the Hamiltonian, using Eq. 10, for a 
given set of transitions in the basis corresponding to 
these transitions such that Htot <^ = E^x {k) (j) , we end 
up with a matrix which is not symmetric, as in C2DEG 
but in addition, here, many matrix elements are com- 
plex as this will be seen in appendix B. To make the 
treatment easier to follow, we adopt the same technique 
as used in Q which consists in writing the hamiltonian 

in a new basis = M <f) where M is a diagonal unitary 
matrix. The new Hamiltonian is then expressed as: 



which is now symmetric and has only real matrix ele- 
ments. 

In the calculations we will neglect the spin splitting A5 
which is small in the case of Graphene 13] but will not 
change anyway the conclusions since the optical transi- 
tions conserve the spin. 

We will furthermore assume in the following that 
there exists some valley splitting Ay first suggested by 
Gusynin et al. [l3| . The existence of such a valley split- 
ting, has been recently supported by different models. 
Some of these models include different types of electron- 
phonon interactions [isl . [l6| : they all predict a linear 
dependence of Ay with the magnetic field. Another one 
|17| invokes EEI with strain induced gauge field yielding 
to a valley splitting which varies like \/B. We assume 
here, for convenience, that Ay is larger than Ag in such 
a way the electrons remain in the same valley ( here the 
valley K for instance) for filling factor v < \. This is 
not necessary true since experimental results [l8| tend to 
favor a situation, where for v <^ \, the system becomes 
spin polarized. The splitting Ay will not be included 
in the present calculations but its consequence will be 
discussed in each case where its contribution could be 
relevant. All the following results for energies are given 
in units of E^ and as a function oi K = I fc I • 



III. MAGNETOPLASMON ENERGIES FOR 

V <2 



One electron optical transitions in Graphene 
(a) (b) 
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(13) 



FIG. 1: (Color on line) Schematic diagram of one-electron 
transitions used in the magnetoplasmon model for 1 < < 2. 
On the left part of the figure (a) are shown the transitions 
implying the n = LL. On the right part of the figure (b) are 
shown the first interband (electron-hole) transitions from the 
n — —2, —1 LL to the n = 1, 2 LL respectively. For < 1, 
the transition a"*" in (a) disappears and the new transition 7" 
appears. The splitting of LL n mimics the spin splitting. 

For z/ < 2, the typical transitions to be analyzed are 
displayed in Fig.l . There are two kinds of one elec- 
tron transitions, in each valley, those which imply the 
n — LL (Fig. la) and those which correspond to inter- 
band (electron- hole) transitions (in Fig. lb are only rep- 
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resented those from LL n = —2,-1 to n = 1,2). As 
already mentioned in the preceding section one can treat 
independently the one-electron transitions implying the 
n = LL and those involving the interband transitions. 

A. Magnetoplasmon energies for transitions 
implying the n — LL 

The different one-electron energy transitions which are 
considered in this case are displayed in Fig. la ior 1 < i/ < 
2. In this figure the splitting of LL n mimics the spin- 
splitting As for clarity but, as already said, this splitting 
is not taken into account in the present calculations. The 
Hamiltonian to be solved is therefore a matrix of rank 5 



written, first, in the basis (j) —{a'^,a ,/?"'", /3 ,7"'") (see 
Fig. la for notations) and then transformed according to 
Eq. 13. The corresponding diagonal matrix M is denoted 
here M0i<,y<2 and has the following diagonal elements: 

M0i<.<2 = {^,e-^^, e-^^ e-'^ (14) 

where Jq is the partial filling factor of the spin-up n = 
LL and the polar angle of the cxciton wave vector. 

For 1 < V < 2, the matrix t)e diagonalised is 

expressed as: 



77° - 




/o'^ioi 

^101 
^101 



/o'(l - /o')(^ooii - -Eooii) y'l - fdVooii \l 



/o K)101 
Vbioi 
/133 
Vbioi 
1 — /o^Vboii 



./o K)101 
Vbioi 
Vbioi 

y'l ^ /o^Vboii 



y^/o^(l - ./o^)(K)oii - -Boon) 
y^l ^ ./o^K)oii 

^ /o^K)011 



/155 



r 



(15) 



The different matrix elements are given in appendices A 
(Eq. Al) and B (Eqs .Bl, B2, B3). 

One can note that the eigenvalues of -ffi<^<2 are identi- 
cal for fQ= or 1 (that is v—\ or 2) whereas those for non 
integer values of u are symmetric with respect to 1/ = 1.5. 

For V <1 the corresponding Hamiltonian ^^o<i><i ^'^ 
be written in the basis (p — {a^ , f3^ , (3^ (see 
Figla), replacing by Jq with a new diagonal matrix 
M denoted now AfOo<iy<i which has the following ele- 
ments: 



MO, 



'0<i/<l 



-lip -I'p I'p 



1-/0" 



=} (16) 



where /q is the partial filling factor of the spin-down 
n = LL. 

The corresponding expressions for the matrix elements 
are given in appendices A (Eq. Al) and B (Eqs .Bl, B2, 

B3). It turns out that the eigen values of H^^^^i are 

symmetric of those obtained for -ffi<^<2 with respect to 
11 = 1. If we adopt a model where Ay is smaller than A5 
for V < 1^ we obtain an Hamitonian with the same eigen 
values which shows that the MP results do not depend 
on this assumption. 

Results for the MP dispersion curves are displayed in 
Fig. 2. For v = 1 or 2 one obtains, for the dispersion 



curves, a solution Ed{K) three times degenerate and one 
solution Eu{K) which have the following analytical ex- 
pressions: 



Ed{K) = + Ci + -ao - Emw{K) 
Eu{K) = Ed{K) -f 41/0101 (if) 



(17) 



where E'lo = Ei—Eq = 2.77xe^ /{kIb) is the one-electron 
energy for these transitions and Ci defined in appendix 
A (Eq. A2) is a quantity a priori divergent which will be 
discussed in section V. 

For non integer values of i^, the solutions Ed{K) re- 
main twice degenerate, the high energy solution remains 
close to Eu{K) and two new solutions appear. The linear 
dispersion near if ~ for Eu{K) is due to the RPA con- 
tribution entering Eq. 10. As compared to the solutions 
found in C2DEG ior v = 2 0, 0] , this contribution is the 
same whereas that of the exciton binding energy is dif- 
ferent. The solutions for K ~0 will be further discussed 
in section V below. 

Following the lines of Ref. 11 we have also calculated, 
in the frame of the MP picture, the optical conductivity 
(see appendix C, Eq. C4) which predicts that Eu{K) 
should be optically active in both polarizations of the 
light (note that the optical vectors are proportional to 

The MP model has been derived without including the 
valley splitting Ay: if such a splitting is introduced, for 
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Graphene: transitions n = -1 ,0 to n = 0,1 



Graphene: transitions n = -2,-1 to n = 1,2 
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FIG. 2: (Color on line) Variation of the magnetoplasmon en- 
ergies in units of /{kIb) as a function of fc/s for the tran- 
sitions involving the n = LL and filling factors < 2. The 
dotted circles denote the degeneracy of the transitions. 



the 71 = LL, we expect a corresponding splitting of the 
optical transition independent on the relative magnitude 
of Av^ and A5. 



Magnetoplasmon energies for transitions from 

the n = -2, -1 to n = 1, 2 LL 



We discuss now the case of interband transitions dis- 
played in Fig. lb. There are, in this case, eight possible 
one-electron transitions and the Hamiltonian is written 
first in the basis (j) ^ {I^, J^, I]^, , J^, ,1^, , J^,) 

and then transformed according to Eq.l3. For < < 2, 
the corresponding diagonal matrix AI denoted here as 
Af-^o<i^<2 has the following diagonal elements: 



(18) 

For these transitions the symmetrized excitonic Hamil- 



12 

<i'<2 



and HIq^^^-^ have matrix elements 



tonians HI^ 

which are given in appendices A (Eqs. A3, A4) and B 
(Eqs. B4, B5, B6). 

The dispersion of MP energies, in units of Ec, are dis- 
played in Fig. 3 as a function of K. The correspond- 
ing one electron energy for these transitions is EI12 — 
El - E_2 = {V2 + l)Ew = 6.69 x 67(^/5). 

As for the preceding case, the solutions are identical 
for z/ = 1 or 2 and symmetric with respect to v = 1. 
For integer values of i^, the eigen-values of the Hamilto- 
nian can be expressed analytically and arranged in two 
groups: (i) two single solutions EI^^ (K) displayed in 
the left part of the Fig. 3 and (ii) two other sets of so- 
lutions EI^^^iK), three times degenerate (see dotted 
circles in Fig. 3), displayed in the right part of the figure. 
They are expressed as: 



FIG. 3: (Color on line) Variation of the magnetoplasmon en- 
ergies in units of /{tds) as a function of Mb for the tran- 
sitions n = —2, —1 to n = 1,2 LL and filling factor v = 1 
or 2. The dotted circles denote the degeneracy of the tran- 
sitions. The magnetoplasmon curves in the left panel are not 
degenerate and are the only optically active transitions. 
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Graphene: transitions n = -3,-2 to n = 2,3 
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FIG. 4: (Color on line) Variation of the magnetoplasmon en- 
ergies in units of 6^/(^/3) as a function of Hb for the tran- 
sitions n = —3, —2 to n = 2, 3 LL and filling factor v = 1 
or 2. The dotted circles denote the degeneracy of the tran- 
sitions. The magnetoplasmon curves in the left panel are not 
degenerate and are the only optically active transitions. 



EI+' (AO =(\/2 + l)(Sio + Ci) + AC2 

+ 4V_12-12 — -E_122-1 



±\l{^? + {'^V-l2-12-E^ll2-2? 



(19) 



and 
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Magnetoplasmon energies for 2 < u < 4 



EI+'-(K) = {V2 + l){Ew + Ci) + AC2 



-E. 



± 



^)2 + (£;_n2-2)2 



(20) 



where all matrix elements entering Eqs. 19 and 20 are 
function of K and given in appendices A (Eq. A4) and 

B (Eqs. B5, B6). Only the solutions EI^'^{K) are 
optically active (see appendix C, Eq. C5). For non- 
integer values of the filling factor the results are very 
close to those presented in Fig. 3 except for two solutions 
of the two groups of degenerate transitions which are no 
longer degenerate for ii' ~ 0. 

In contrast to the case of transitions implying the n = 
LL, a splitting of the transitions equal to ao/S and due 
to electron-electron interactions is expected for i^T ~ 0. 

Note however, here, that the introduction of a valley 
splitting Ay should only provide an additional compo- 
nent either linear in B or in ^fB depending on the origin 
of this valley splitting. 



Magnetoplasmon energies for transitions from 

the n = —3, —2 to n = 2, 3 LL 



In this case the corresponding Hamiltonian HI^^^^2 

has the same structure that HIq^^^2 ^-nd therefore only 
the values of matrix elements are different. They are 
given in appendices A (Eqs. A5, A6) and B (Eqs. B7, 
B8). The dispersion of MP energies, in units of E^., 
is displayed in Fig. 4 as a function of Mb- The cor- 
responding one electron energy for these transitions is 
EI23 = E2- E_32 = (V3 + V2)Eio = 8.72 x eV(K/s). 
The solutions are formally identical to those given in Eqs. 
19 and 20 with the appropriate changes for the matrix el- 
ements given in the appendices A (Eq. A5) and B (Eqs. 
B7, B8). The splitting of the transitions for iiT ~ is 
here equal to ao/16. 

We, now, evaluate the exciton energies for > 2. 



IV. MAGNETOPLASMON ENERGIES FOR 

2 <u <6 



In this case we have to treat the problem depicted in 
Fig. 5 for the one electron energy transitions. Note that 
we have here two types of transitions those implying the 
n = LL and those between LL n = 1 and n = 2. 
Because the corresponding one electron energies are dif- 
ferent they are treated independently. 



One electron optical transitions in Graphene: 2 <v< 4 
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FIG. 5: (Color on line) Schematic diagram of one-electron 
transitions used in the magnetoplasmon model for 2 < ;/ < 4. 



For 2 < ;/ < 3, we have to write, first, the Hamiltonian 
in the basis 4> — {a^ , , /3+ , , 7^ } and for 3 < < 4 

in the basis (f) = a^, 7^, 7^} and transformed 
them according to Eq. 13. For 2 < < 3 and 3 < 
< 4, the corresponding diagonal matrices M, denoted 
here as M122<i/<3 and M123<y<4 respectively, have the 
following diagonal elements: 



M12 



2<iy<3 



-{1,1,1, 



l-/r \ f : 



M12 



3<iy<4 



{1,1, 



1 - /i+ \ .ft 



-} (21) 



We will concentrate the report for filling factors 2 < 
v < 6. The contributions of exchange are given in ap- 
pendix A (Eq. A7). It turns out that the problem to 
solve is symmetric with respect to v — 4 and therefore 
we will detail the treatment for 2 < < 4 and will note 
only the main changes for A < v < 6. 



where and f^^ are the partial filling factors of the 
spin-down and spin-up n — 1 LL respectively, attached 
to the valley K' with our convention. One gets for the 

corresponding Hamiltonians -ff2<iy<3 -ff3<y<4 the fol- 
lowing expressions: 
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and 



rrl2 

-"2<;y<3 



TT12 

-"3<i^<4 



M]101 

1 ^ /i Vbioi 




hll 
Vbioi 
1 — /I'Vbioi 





Vbioi 



/i Vbioi 




VoiOl 



1 — /i^Vbioi 





K)101 

Vbioi 

/i33 



1 ^ /i Vbioi 




1 ^ /i^Vbioi 

1 ^ /i^^lOl 
^33 






1 - fiV^ - 







1 - /fVbioi 
1 - fiVnoi 



^44 









/l44 
/l^^l212 





^55 







/l^Vi212 
^55 



(22) 



(23) 



where the new matrix elements entering these matrices 
are given in appendices A (Eqs. A7) and B (Eqs. B9, 
BIO). The resuhing dispersion curves are displayed in 
Fig. 6. The corresponding one-electron energies for both 
types of transitions are i?io — 2.77 x /{k,Ib) for the 
higher ones and £'12 — E2 — Ei — 1.15 x e^/ {kIb) for the 
lower ones. 
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FIG. 6: (Color on line) Variation of the magnetoplasmon en- 
ergies in units of /{kIb) as a function of Mb for different 
filling factors 2 < v < 4. The dotted circles denote the degen- 
eracy of the transitions. 



The eigen-values of the Hamiltonians (Eqs. 22 and 
23) can be expressed analytically. For 2 < v < ?> there 
are two solutions identical to Ed{K) (Eq.l7) which re- 
main degenerate and which are not optical active. Two 
other non degenerate solutions denoted as EQ^^^^z^K) 
are given by: 



EQt<u<z{K) = -[hll + hi4 + 2Voini 



(24) 



±J{hM - hll - 2^0101)2 + 12yo'ioi(l - /r)] 



and a third one E122<ciy<3{K) — hc,c,[K). These transi- 
tions are all optically active with a relative weight de- 
pending on the filling factor. 

For 3 < < 4, it remains one eigen- value solution 
identical to Ed{K) (Eq.l7) and two groups of optically 
active non degenerate solutions denoted as EO^^^^^{K) 
and E1^^^^^{K) with the following analytical expres- 
sions: 



EO 



3<u<4 



(K) 



[hll 



'33 



2Vn 



0101 



±\J{h33 - hll - 2Foioi)2 -H8FoV(l - /i+)] 



(25) 



and: 



El 



3<iy<4 " 2^h44 



'55 



(26) 



±V (/144 - h5 



4V^1212/i+ 



for which all matrix elements arc function of K. 

The corresponding optical vectors for these transitions 
are given in appendix C (Eqs. C7 and C8). 

B. Magnetoplasmon energies for 4 < 1/ < 6 

It is easy to see that, in this case, the structures of 
the corresponding Hamiltonians H^^^^^ and H^^^^g are 
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FIG. 7: (Color on line) Variation of the magnetoplasmon en- 
ergies in units of /{k,Ib) as a function of kls for different 
filling factors 4 < < 6. The dotted circles denote the de- 
generacy of the transitions. 



We will not discuss the case of interband transitions for 
this configuration of filling factors but their correspond- 
ing Hamiltonians HI^'^j^^q and H^^^^q are modified in 
two respects: for both of them the exchange contribu- 
tions entering the diagonal elements are different and for 

^^2<iy<6 transitions are now filling factor dependent 
in such a way the corresponding transitions disappear at 
v = 6. 

The case of filling factors 6 < < 10 and following 
ones will not be discussed as well but the corresponding 
treatment is formally similar to the case 2 < v < 6 with 
a different set of one electron energies, exchange contri- 
butions and matrix elements. 

We are now focussing the discussion on the results ob- 
tained for ii' ~ which could be compared to magneto- 
optical absorption measurements. 



DISCUSSION OF THE RESULTS FOR ~ 



symmetric with respect to those given in Eqs. 22 and 23. 
At present, this is the n = 1 LL attached to the valley 
K (in our convention) which starts to be filled and the 
notation refer to this LL. Of course some of the diago- 
nal matrix elements are changed but results are formally 
similar and the corresponding exciton dispersion curves 
are displayed in Fig. 6. All the eigen- value solutions of 
Fig. 6 can be expressed analytically: 

For 4 < 1/ < 5 one gets one solution Ed2{K) — h^^{K) 
which is not optically active and two groups of optically 
active solutions: 



For K = I fc Z_b| — 0, all the Hamiltonians are reduced 
to their diagonal elements which are given in appendix 
B. The reason is that all off-diagonal elements are pro- 
portional to K or . We will call the corresponding 
solutions, at i<r ~ 0, E^j^^ and £'j|"p"^^' for intra-LL 
transitions and interband transitions respectively. We re- 
strict the discussion to those solutions which are optically 
active. All results are function of ao = ■7^^/^ = 0.627 in 
Coulomb units. We then get the following results: 

For the transitions i?°)p which involve the n = LL: 



'-22 



(27) 



and: 



El 



4<i/<5 



:[h3 



Vl21 



(28) 



± J{h4i - h33 + 1^1212)2 + SV^^J^] 



for which all matrix elements, dependent on K, are given 
in appendix B (Eq. Bll). 

For 5 < J/ < 6 one gets two solutions Ed2{K) = 
h44{K) (same expression as for 4 < < 5) which 
are not optically active, one optically active solution 
E05<ciy<e{K) = /ill (if) and two other optically active 
solutions denoted as EI^^^^^q. 



El 



± 

5<i^<6 



*33 



2Vu 



(29) 



±J{h 



^33 



2Fl212)2 + 12V,1,J + 



The related matrix elements, dependent on K, are given 
in appendix B (Eq. B12). 



Ci 



(for V = 1,2,3,4,5) 



i?io + Ci + ^(-3 + ^(6-j.)) 



(for 5 < J/ < 6) 



(30) 



As clearly apparent in Fig. 6 and 7, this transition is 
split for non integer values of > 2. The high energy 
component of this split level has an energy which in- 
creases with v but its oscillator strength decreases with 
V going to zero at integer value of v. 

For the transitions E^ip which involve the transitions 
between the n = 1 and n = 2 LL: 



Efip = (V2-l)(^io4 
-^(1 + 2V2) 
(fori/ = 3, 4, 5, 6) 

E]^,P = (V2-l)(^io4 

(for 2 < i/ < 3) 



Ci) + AC2 



Ci) + AC2 
-2V2){i^-2)) 



(31) 
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In this case also this transition is spht for non integer 
values of v. The high energy component of this split level 
has an energy which increases with v but its oscillator 
strength decreases with v going to zero at integer value 
of V. 

The optical active interband transitions Elj^p which 
involve the transitions between the n = —2,-1 and n — 
1,2 Landau levels are split by an amount ^ but the 
mean energy Elmfjp = iEI^{0) + EI^{0))/2 (Eq. 19) 
has the following expression: 

Elmf.p = (V2 + l)(^io + Ci) + AC2 - — (32) 

whereas the mean energy for the optically active transi- 
tion implying the n = —3,-2 and n = 2, 3 Landau levels 
is expressed as: 



TABLE I: Evaluation of the re-normalized velocities c, at in- 
teger values of the filling factor, for difi^erent transitions. 



Transition 
n to TTl 


V 


C t^iU TTl/ S ) 


C l^iU TTl/ S) 


o to 4 


0, iU 




u.yy ± U.Uz 


2 to 3 


6,8 




1.01 ±0.02 


1 to 2 


4,6 




1.04 ±0.02 


to 1 


2,4 


1T2±0.02 " 


1.12 ±0.02 




2 


1T2± 0.004 ' 




-1 to 2 


2 


1.18 ±0.02 " 


1.16 ±0.02 


-2 to 3 


2 




1.16 ±0.02 


-3 to 4 


2 




1.16 ±0.02 



"from Ref. |13| 
''from Ref. 



Elm'j^P = (x/3 ± V2)(£;io + Ci) + AC3 - —ao (33) 

256 

For a given value of the Fermi velocity vp ( here equal 
to 0.86 X lO^m/s) the energy £^10 is determined (here 
EiQ = 2.77) and in Eqs. 30 to 33 the only unknown 
parameter is Ci. 

As already said this quantity defined in appendix A 
(Eq. A2) is divergent. The occurrence of such a prob- 
lem is not specific of the Graphene properties because 
it is also present in C2DEG though it was not explic- 
itly formulated. The reason why this term appears here 
is that we wanted to define a common energy scale 
for intraband and interband transitions. This will cor- 
responds in C2DEG to impose a common energy scale to 
cyclotron-type transitions and inter-band excitonic tran- 
sitions. There was an attempt to treat this later tran- 
sitions in GaAs, in an another context, but using the 
same theoretical model [2l[ and indeed the same problem 
of divergence of the exchange interaction among the va- 
lence band levels was found without being able to solve it. 
Therefore this problem is not specific to Graphene but, 
in that case, one can solve it, at least, in a semi-empirical 
way. 

The divergence of Ci is due to the infinite summa- 
tion over LL (see Eq. A2) which is physically artifi- 
cial. We could then define, as was done in Ref. ^8^, 
a cut-off value on energy or number of LL, but this 
limit is quite arbitrary. We propose to treat the prob- 
lem in a semi-empirical way, using Ci as a parameter 
fitted, for one type of transitions, to experimental data 
and then to deduce all the re-normalized velocities at- 
tached to the other transitions. Doing so, we implic- 
itly assume that all the re-normalization of the veloc- 
ity, for the fitted transition, is only due to electron- 
electron interactions neglecting any possible contribution 
from electron-phonon interaction which may be impor- 
tant in carbon based compounds. Among experimen- 
tal data which could be used for this fitting, those re- 
lated to magneto-transmission measurements [ill . [l2j are 



those which are expected to reflect the magneto-plasmon 
picture developed is this study. Another set of data, 
based on photoconductivity measurements on ex-foliated 
Graphene [13] can also be considered to compare re- 
sults. We will use the data of Ref.[ll|, obtained on 
ex-foliated Graphene, to fit Ci to the E'^p transition 
at V = 2. In this reference, the re-normalized Fermi 
velocity cgf = (1.12 ± 0.02) x lO^m/sec which, from 
Eq.30, leads to a value Ci = 1.31 ± 0.06. When in- 
jecting this value in Eq. 31, we predict a re-normalized 
Fermi velocity c/J^ = (1.163 ± 0.02) x lO^m/sec to be 
compared with the corresponding experimental value (l^ 
c/e^ = (1.18 ± 0.02) X 10*^TO/sec. The agreement is rea- 
sonable. As seen on both experimental and theoretical 
grounds, the re- normalized velocity differs for different 
transitions. One can then try to evaluate these veloci- 
ties for other transitions. The results are given in Table 
1 for integer values of v. The value quoted in the ta- 
ble, from Ref. [23| . corresponds to the transition n = to 
n = 1 whereas the corresponding value for the transition 
n = -1 to n = is found to be (1.07±0.004) x lO^m/sec. 
This corresponds to an asymmetry of the conduction and 
valence levels not taken into account in our model but 
also not reported in Ref. [13] for this transition. 

In table 1, we have included the results of the re- 
normalized Fermi velocity for transitions between LL 
n = 3ton = 4, n = 2ton = 3, n = —3 to n = 4 
for which we have calculated the diagonal elements of 
the corresponding Hamiltonians. 

As clearly apparent, from Table 1, c varies strongly 
with the transition though it seems to be relatively con- 
stant for all interband transitions. This qualitative fea- 
ture is also observed in experiments performed on epi- 
taxial multi-layer Graphene [ll|, |2^. However in these 
experiments the reported values of c"^^ for all transitions 
is the same and equal to (1.03± 0.01) x lO^m/s which is, 
at present, not understood. 

One could, a priori, think that the results obtained 
are dependent on the value of vp adopted in the calcula- 
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tions. In fact one can vary this value over a large range, 
for instance, from 0.80 x lO^m/s to 0.90 x lO^m/s getting 
values for Ci = 1.50 to 1.18 respectively but the quan- 
tity which enters the Hamiltonians is in fact Eqi + Ci 
which remains constant, independent on vp and equal 
to 4.08. This value has been adopted to calculate the dis- 
persion curves of Figs. 2, 3, 4, 6 and 7. It is therefore not 
possible from experimental data on energies to determine 
Ci but all transitions are now given with a common en- 
ergy scale. On the other hand the oscillator strengths 
of the transitions are proportional to Vp and then abso- 
lute transmission measurements could in principle give 
information on vp and therefore on Ci. 



VI. CONCLUSIONS 

In conclusion, we have developed . within the Hartrec- 
Fock approximation, a full treatment of the magnetoplas- 
mon picture in Graphene valid for a very large range of 
magnetic fields. This model, applied for filling factors up 
to 6, shows that the electron-electron interactions induce 
different effects: (i) for some of the transitions these in- 
teractions lead to a splitting of the optical transitions and 
(ii) they are responsible for a strong re-normalization of 
the Fermi velocity as observed in magneto-optical exper- 
iments. This re-normalization is found to be dependent 
on the type of investigated transitions. The optical con- 
ductivity components have been evaluated showing that 
the oscillator strength of the optical transitions is propor- 
tional to Vp and not to the square of the re-normalized 
velocity. The theory has been derived for all transitions 
with a common energy scale which should allow a direct 
comparison of its predictions with future experimental 
works. 



1. Exchange contributions to H^^^^2 ^"^d J?o<i'<i 
For -ffi<y<2 obtain: 

^^(-E'O,m,O,Tn(0) — -E'l,m,l,m(0))/m = 
m 

^ao(2/o+-l) + Ci 

^(-E^O,m,0,m(0) — Ei^m,l,m{0))fm = J^fO + C*! 
m 

■ 

-Jao(2/+-l) + Ci (Al) 

where: 

V2^7o Vm-t-1 

and are Laguerre polynomials of argument x'^ in 
this Appendix. 

The quantity Ci diverges due to the simplifying as- 
sumption of the infinite linear dispersion of the Graphene 
bands. The summation has to be truncated at some level 
or this parameter has to be fitted to experimental data 
(see section 5). 

For i?o<i/<i obtain for the exchange part the same 
expressions than those given in Eq.Al replacing /q" by 
fo- 
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2. Exchange contributions to J?i<„<2 ^( 
For HIl^^^2 0116 gets: 



APPENDIX A: EXCHANGE CONTRIBUTIONS 

We report in this appendix the explicit expressions for 
the contribution of the exchange energies entering the 
diagonal elements of the different Hamiltonian matrices 
in units of Coulomb energies. We introduce the notation 
ao = ^ \/¥ wliich characterize the exchange interaction 
in C2DEG at = 1. Applying the expression given in Eq. 
12 we obtain successively the contribution of exchange 
for the different Hamiltonians. To simplify the notations 
we will drop the superscript ~ from En,m,n,m{0) meaning 
that all these quantities are real. 



m 

^(-£^-2,m,-2,m(0) — -E^l,m,l,m (0))/^ = + ^^12 

m 

Y,iE-l,m,-l,m{0) - ^2,m,2,m(0))/+ = 

m 

^(£_2,„,-2,m(0) - Ei,rnA,mmfiL = 
m 

-^(2/o+-l) + C/i2 (A3) 
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where: 



-y I 

* TO=0 " 



2 r2 

d.:e-^^^L^[l + :%] 
Vm + 1 \/2 



= (\/2 + l)Ci+AC2 



- 00 



m=0 



2 a; 



a/to + 1 



(A4) 



C/12 in this equation also diverges like Ci but AC2 con- 
verges to a value -0.156. 

Similar expressions hold for -ff-fo<i^<i when replacing 



3. Exchange contributions to ifi<^<2 and ffo<i'<i 



For HIf^^^2 oiie gets: 

^(-E-2,m,-2,m(0) — -B3,m,3,m(0))/r 

m 

^(-E'-3,m,-3,m(0) — -E2,m,2,m (0))/~ 
m 

^(^_2,™,-2,m(0) - £;3,„,3,m(0))/+ 
m 

^(2/o+-l) + C/23 

3,m, — 3,m 

(0) - E2 ,m,2,m (0))/+ 

m 



32 



23 



32 



23 



32 



(2/0+ - 1) + CI: 



23 



where: 
C/23 

AC. 



1 /■'^ 



dxe ^ 



^2 r 1 r 1 

r 1 I -^2 N 

"^x/2 ^/3^ 



= (\/3 + \^)Ci +AC3 



^ 00 

4i: 



m=0 



\/m+ 1 



(A6) 



C/23 in this equation also diverges like Ci but AC3 con- 
verges to a value -0.467. 

Similar expressions hold for HIq^^^-^ when replacing 
/o+by/o". 

Comparing Eqs. A4 and A6 one can formally extend 
the treatment and find that, for any interband transition 
from LL —ptoLL q=p+l the corresponding divergent 
term CIpq entering the exchange contributions is given 
by CIpg = {y/p + y/q)Ci) + Fp^g where Fp^g is finite. 



4. Exchange contributions to i?2<i/<6 

The different contributions to the exchange for the dif- 
ferent Hamiltonians are: 



^(■E^O,m,0,m(0) — El,m,l,m{0))fm 
■ 

«o(^-^/?) + Ci 

^(-£'l,m,l,m(0) — -£^2,m,2,m(0))/m ~ 
m 

g(5 + (26-4^/2)/±)+Cl2 (A7) 

where C12 = (^2 - l)Ci + AC2. 

Here also formally, when extending the treatment, one 
finds that, for any intra-LL transition from LL p to 
LL q = p + 1 the corresponding divergent term Cpg 
entering the exchange contributions is given by Cpg = 
^ + Gp^q where Gp,g is finite. Note however 

that in general Gp^g is different from Fp^g. 



APPENDIX B: HAMILTONIAN MATRIX 
ELEMENTS 

We report in this appendix the explicit expressions 
for the matrix elements of the Hamiltonian matrices in 
units of Coulomb energies. Results arc given as a func- 
tion of K = \klB\- For simplicity we adopt the same 
notation hij for noting the matrix elements of all ma- 
trices but their expression is specific of the case under 

consideration. All matrix elements V 



n-i ,n2 ,n3 ,714 



ill) and 



(A5) E„ ( /c ) are evaluated using Eq. 11. 



1. Matrix elements of H^^^^^ and ^fo<i/<i 



For i?i<y<2 obtain: 



/in = + ^ao(2/o+ - 1) + /o+(T^oioi - -Eoiio) + Ci 

/l22 = /l33 = /l44 = £^10 + ^ClO + (K)101 — Equq) + Ci 

/i55 = - \ao{2f+ - 1) + (1 - /o+) 



X {V- 



10-10 



- E. 



100-1 



) + Ci 



(Bl) 



The matrix elements V„i„2„3„4 entering the Hamilto- 
nian i/i<y<2 



(B2) 
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with Vo, 1,0,1 = Vqwi and Vo, 0,1,1 = Vboiie^*"^ where is 
the polar angle of the exciton wave vector. 



The matrix elements En 
Bl are: 



entering Eqs. 15 and 



2[*(2'l;-^)-4*(2'l;-^)] (B3) 



Eomi{K) = — 



3^ 
~32" 



2^(2'^'-^) 



where -Eq^i.i.o = -Bono, -^0,0,1,1 = -Booiie^"^ and 
$(a, 6; z) is the confluent hypergeometric function. 

For -ffo<i^<i matrix elements arc identical to those 
given in Eqs. Bl, B2 and B3 replacing by /J" when 
appropriate. 



2. Matrix elements of HIl'^^^^ ^iid ^?^o<i/<i 



The matrix elements Enin2n3n4 entering B4 are: 



/¥ 1 K"^ 

-$Ai;-^) + V;i;--)--<J>(-;i;--)] 

^2 2 ^ 4 ^2 2 ' 64 ^2 2 ^' 

p i^^V^. 3(l + x/2)^ 5.,. 

/7^x;^/7 „ 105\/2^,9 „ 

- 15(2 + ^/2)$(-; 3; --) + -^H^, 3;--)] 



(B6) 



where £'-1,2,2,-1 = -B-122-1 and £'-1,1,2,-2 = 

£-112-26^*'''. 

For HIq'^jj^^ two columns of the matrix Hll^^,^2 are 
inverted but the eigen values are the same with fQ re- 
placing /o" in Eq. B6. 



The matrix elements hij of HI\'^^^2 are: 

hii = EI12 + 7^+^2-1- V-12-12 — £-122-1 
16 

/l22 = £/l2 — 777 + ^"2 + ^-12-12 — £-122-1 

16 

^33 — EI12 + (^■^0'" — 1) + C2 + ^-12-12 — £-122-1 
/l44 = £-^12 — '^i'^^Q — 1) + C2 + V-12-12 — £-122-1 
/l55 = £/l2 + 7-^ + C2 + V-12-12 — £-122-1 

16 

/l66 = £^2 — TT + (^2 + V-\2-l2 — £-122-1 

16 

/l77 = £/i2 + — + C2 + ^-12-12 — £-122-1 

16 

^88 — EI12 — — + C2 + V-\2-\2 — £-122-1 

16 

^12 = ^34 = ^56 = ^78 = ^^-11-22 — £-112-2 
/il3 = ^15 = ^17 = ^24 = /l26 = ^28 = ^-12-12 
^35 = /I37 = /l46 = /l48 = /l58 = /i68 = 1^-12-12 
/ll4 = /I16 = h\% = /l23 = /l25 = /l27 = V-W-Tl 

hs6 = /i38 = /145 = hir = /157 = her = V-n-22 (B4) 

where £/i2 = (v^ + l)£io- The matrix elements 
Kinsnan* entering Eq. B4 are: 



y-i2-i2(i^) =y-ii-22(i^) = ^e-'^[(3-2y2) 
+ {V2-2)- + -] 



(B5) 



with y-1,2,-1,2 = ^-12-12 and V^-i, 1,-2,2 = 
F-ii-22e2*^. 



3. Matrix elements of HI^^^^^ ^"d H^^^^^ 



The matrix elements hij of i?/i<y<2 are similar to 
those given in Eq. B4 when replacing V^-12-12, V-n-22, 

£-122-1, £-112-2 by V-23-23, ^-22-33, £-233-2, 

£-223-3 respectively and £/i2 by £/23 = (•\/3 + ^/2)£lo• 
The new matrix elements En-^n^nsm are here: 



£_ 



233-2 



{K) = 



15^,5 ^ K^ , 795^,7 ^ 

H $(-;l; ) $(-; 1; ) 

4 4'' 2 ' 192 4' ' 2 ^ 



315^,9 , 
+ 128^^2^^^ 



^ 1536^ 2 ' ' 2 



£-223-3 W = -^^[3(5 + 2V6)a>(| 3; 

-15(4+^).4;3;4) 

.105(^^).(^;3;4) 



2 ^ 



(B7) 



where £-2,3,3,-2 = £-233-2 and £-2,2,3,-3 = 

-^'-223-36 ^• 

and the corresponding matrix elements Vmn^nim'- 



13 



,3 7 



(^-6) , ' +(l-/r-)(Vo.o.-W + C, 

^ ^ 2' ^i2+^gU+ 2 (B12) 

with V^_2,3-2,3 = V^-23-23 and V^-2,2-3,3 = + /i+(Vi2i2 - £^1221) + 

^-22-336^''^. 

^33 = - /l55 - £^12 + ^(21 - 10 V2) 

4. Matrix elements Of + (1^212 - ^1221) + 

with the corresponding new matrix elements entering Eq. 

For Hll^^^ we obtain: B9, BIO, Bll and B12: 

hii = h22 = /»33 - ^10 + ^ao + (M)ioi - Soiio) + Ci 1^212 W = "g-e-— [1 + V2 - 

h,, E,o + ao{^ + (1 - /r)(^oioi - i?oiio) + Ci = i?-i22-i W (B13) 
4 o 

/^55 = i?12 + ^(1 + ^^'f^ )/r ^l^^'^^ ^1.2.1.2 = V^1212 

16 6 



+ /l (^1212 ~ -E1221) + C2 



(B9) 



For ffi2,<4 : 



APPENDIX C: OPTICAL CONDUCTIVITY 

In this appendix we derive, following the lines of 
g Ref. 0, the corresponding expressions which allow to 

hii — = Eio H — ctQ + (Vbioi ~ -£'Oiio) + C*! calculate the optical matrix elements of the MP curves 

^ which enter in the optical conductivity a[hhj) which has 

/133 - i^io + "0(7 - Ift) ^^"^ components: 

/i44-^12 + ^(21-10a/2) (BIO) Lo ^ {Ei^pf - {hLof 

+ {V1212 - E1221) + C'^ ^ ^GB ^ '2rMMl-L'i}{Ml-L'jY 

_ a, 19-10^/2 " * ^ ^ {Elpf - {h^Y 

n^b - E12 + —(1 ^ j/i 

+ /i^(^i2i2 ^ £^1221) + C2 where the summation is performed on all MP transitions 

of energy -B^^p with the corresponding eigen vector Lj . 
For -ffl<^<5 : In Eq. CI, Gb 1/{2t:{IbY) is the density of states 

of a single LL, A^n and M± are optical vectors with 

/ill = i?io + 7ao + (^0101 - £^0110) + Ci components M-^^T'l^^ and AC/n;a, respectively. M 

4 is the matrix used to symmetrize the Hamiltonian and 

/122 = Eio + q;o(- — -/ ) depends on the set of transitions which are considered 

4 8 ^ (see Eqs. 16, 18, 21). aj denotes one of the transition 

+ (1 — /i^)(Vbioi — -Bono) + Ci n to m belonging to this set of transitions. J-a;aj with 

Q,g X9 — lOv^ (Bll) ^ II defined in reduced units {K standing for k Ib 

ft-33 = Bi2 + — (IH ^ )/f and u for a;/^^) as: 

+ A (^1212 — B1221) + C2 

^44 = /^55 = Bi2 + g(21-10V2) .F„;™,„(]?) = j due^^^-[{F^(u+^)rnF:,{u-^))] 

+ (^1212 - £'1221 )+ C2 (C2) 
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where the function [i^*F_F] denotes the scalar product 
and the velocity operators are: 



K 



VF 







=VF 



e-»^ 




(C3) 



andV;^ = (Kf' )*. 

In the one electron picture the selection rules for opti- 
cal transitions between LL m ^ n are (^iml.lnlii [ill- fl^- 



3. Optical vectors for HI^^^^2 and Hlg^^^^ 

In this case we get for the components of the corre- 
sponding optical vectors M\\ and M±: 



{1,-1,1,-1,1,-1,1,-1} 



,1 I .K^ 



{1,1,1,1,1,1,1,1} 



(C6) 



1. Optical vectors for ^^i<y<2 and H^^^^i 

In the case of H^^^^2 obtain the following compo- 
nents of A^ii and M._l: 



^e-^{V/oM,i,i,-Vi-/o+} 



M\ = ^^e-'^{^f+, 1, 1, -1, - Vl - /o+} (C4) 

whereas ioi < v < 1, has to be replaced by Jq. 

It can be easily verified that for X ~ where all tran- 
sitions become degenerate, there is a sum rule such that 
• £jp — 2vp independent of the filhng factor 

whereas ' Lj){M± ■ Lj)* ~ ivvp. One therefore 

recover the selection rules obtained for the one-electron 
model. Note that the Fermi velocity entering in the op- 
tical matrix elements is that existing in the absence of 
electron-electron interactions. 



Optical vectors for HI\'^^^2 and HI^^^^^ 



Following the same approach we get for the compo- 
nents of the corresponding optical vectors M\\ and M±: 



{1,-1,1,-1,1,-1,1,-1} 



MJ\ 



i—e 4 (l-X^ — ) 

2 ^ 2V2 



(C5) 



{1,1,1,1,1,1,1,1} 



which are no longer dependent of the filling factor for 
i> < 2. It can be shown that the only optical active tran- 
sition are those corresponding to the solutions EI^^ (K) 
ofEq. 18. 



4. Optical vectors for H^'^^^^^ 

Here one gets for the components of the corresponding 
optical vectors M\\ and the following relations where 
we have defined the functions p\\ (K) = --i=(l — "^^"^^^) 
and,,(i.).^(l-i^): 



For 2 < z/ < 3: 



A^||^2<i/<3 



vpe' 



V2 



{1,1,1, ^Jl-fo,P\\iK),/U} 



■Ml 2<^<3 = « 



(C7) 



V2 



-e 4 



{1,1,1, ^Jl-fo,P±iK)y/u} 



For 3 < ly < A: 



vpe 



3<i/<4 



V2 



-e 4 



(C8) 



%/2 



-e 4 



{1, 1, ^1- f+,p^{K),p^{Ky^} 

For A<v <b: 



-lip 



4<i^<5 



x/2 



-e * 



{1, ^l-f^,V\\ [K] Pii [K) , {K)} 



V2 



{1, ^1 - /„- , Pi (X) y^. Pi (if) , Pi [K] } 

(C9) 
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For 5 < 1/ < 6: 



5<u<6 



V2 



p»{K)Jf+,p»iK),p»{K),p»iK)} 



vpe 



V2 

{Jl-f+,p^{K)df, 



p^iK),p^{K),p^iK)} 
(CIO) 

It can be verified that, for X ~ 0, there is a sum rule 
such that for the transitions n = to n = 1, ^ ■ ■ 



whereas for the transitions n = 1 to n = 2, 



We also get EjiM ■l'jKmI- ^)* = i \M\\ ■ h?- 
One therefore recovers the selection rules obtained for the 
one-electron model. 

To conclude this part, it is worth comparing these 
results with those obtained in C2DEG where the in- 
troduction of electron-electron interactions change the 
selection rules 0, Q with respect to those obtained in the 
one-electron picture. In the case of Graphene, the selec- 
tion rules are in general, similar to the results obtained 
in the one-electron picture, except that the oscillator 
strength is condensed into one or two branches of the 
MP curves and that the strength remains proportional 
to v^p and not to the square of the re-normalized velocity. 
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